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ABSTRACT 

Young massive clusters (YMCs) with stellar masses of 10^ — 10^ Mq and core stellar densities of 
10"* — 10^ stars per cubic pc are thought to be the 'missing link' between open clusters and extreme 
extragalactic super star clusters and globular clusters. As such, studying the initial conditions of 
YMCs offers an opportunity to test cluster formation models across the full cluster mass range. 
GO. 253 + 0.016 is an excellent candidate YMC progenitor. We make use of existing multi-wavelength 
data including recently available far-IR continuum (Herschel/Hi-GAL) and mm spectral line (HOPS 
and MALT90) data and present new, deep, multiple-filter, near-IR (VLT/NACO) observations to 
study G0.253 + 0.016. These data show G0.253 + 0.016 is a high mass (1.3 x 10^ Mq), low temperature 
(Tdust ^20 K), high volume and column density (n 8 x lO^cm"'^; N//^ ^ 4 x lO^'^cm^^) molecular 
clump which is close to virial equilibrium (Mjust ~ -^^viriai) so is likely to be gravitationally-bound. 
It is almost devoid of star formation and, thus, has exactly the properties expected for the initial 
conditions of a clump that may form an Arches-like massive cluster. We compare the properties 
of GO. 253 -I- 0.016 to typical Galactic cluster-forming molecular clumps and find it is extreme, and 
possibly unique in the Galaxy. This uniqueness makes detailed studies of GO. 253 -I- 0.016 extremely 
important for testing massive cluster formation models. 

Subject headings: ISM: clouds — ISM: molecules — stars: formation — galaxies: star clusters 



1. INTRODUCTION 

One of the ground breaking discoveries of the Hub- 
ble Space Telescope was that massive stellar clusters, 
with properties that rival (or even exceed) those found 
in globular clusters in terms of mass and st ellar den- 
sity, are st ill forming in the universe today ( jHoltzman 
et al. 1992). These young massive clusters (YMCs) can 
be broadly defined as stellar clusters more massive than 
10'* M© and with ages that are <100Myr but that also 
exceed the current crossing time by a factor of a few, so 
are gravitationally-bound (s ee e.g. Ashman fc Zepf|2001 
Portegies Zwart et al.|2010 ) . The discovery ot YMCs in 
the Galaxy has allowed a detailed study of their struc- 
tural and stellar properties (e.g. Arches, Quintuplet, 
Westerlund 1 and 2, RS CGl, 2 and 3, G LIMPS E-COl, 
NGC 3603: [Figer et al.L19 99; C lark et al.v2 005; Ascenso 
eTaT] [20071 IFiger et ai.|l2006| (Davies et al.||20ll[ |Bec-' 



cari et al. 20101. Like their extragalactic counterparts, 
these Galactic YMCs have stellar masses of 10"* — 10^ Mq 
and core stellar densities of 10^ — 10^ stars per cubic 
pc. The stellar masses and densities are orders of magni- 
tude larger than typical open clusters and comparable 
to those in globular clusters, super star clusters seen 
in merging/starburst galaxies and those inferred from 
observations of massive clusters forming at the epoch 
of peak star formation density (z ~ 2 — 3) (jW hitmore] 
2000l|Johnson fc Kobulnicky|2003[ JSwinbank et al.,201(]f 
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Danielson et al. 2011). As such, YMCs are potentially 
local- universe-analogs of these massive extragalactic and 
high-z stellar clusters. Within our own Galaxy, obser- 
vations show a smooth continuum in the properties of 
young clusters, suggesting that YMCs merely represent 
the extreme end of the open cluster (OC) distribution 
function. Thus, understanding how YMCs form may pro- 
vide a direct handle on i) globular and super star cluster 
formation, ii) cluster formation at the epoch of peak star 
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formation density (2: ~ 2 — 3), iii) the most extreme con- 
ditions for star formation in the local universe, and iv) 
whether or not there are different formation modes across 
the cluster mass range. 

Despite their importance, the formation process of 
YMCs is not well constrained observationally due to the 
lack of known potential gas clouds as progenitors. As a 
result it has not been possible to study the initial condi- 
tion s of the molecular gas from which YMCs form (see 
e.g. [Santangelo et al.|[2009l ). Pro gress requires first find- 
ing and then deriving the properties of the most massive 
molecular gas clouds which are on the verge of form- 
ing YMCs. Such YMC progenitor clouds are expected 
to have a high mass (^IO^Mq) and small radii (a few 
pc) with high volume and column densities (^10'* cm~^; 
>10^'^ cm^^). Before the onset of star formation, the gas 
is also expected to be very cold (~ 10 — 20K). In order to 
resolve the internal structure of objects with radii ^pc 
requires limiting the search to gas clouds within our own 
Galaxy. 

Compared to the objects found in recent surveys 
of dense cluster-forming molecular clumps, o ne object, 
GO. 25 3 + 0.016 , stands out as extreme (e.g. Sridharan 
al][2 002 Hil l et al.||2005l|Beuther et al.|20u!>l IScEuITer 
eTaT 2009 ; Bally et al.l2010''Aguirre et al.'boilVMolinari 
et al.|2 010a; Rathbornc ct al. 2006; Simon ct al. 2006a b; 
Pillai et al.||2006nPeretto fc Fu ller 2009; Jackson etaT 
20061 IRath borne et al.||2009| ILongmore ct al. 2007! IPur- 
cell et al.|| 2006, 2009; Walsh ct al. 2008, 2011; Jackson 
et al.|2011| ). Identified as an infrared dark cloud (IRDC), 
it is clearly seen as a prominent extinction feature in mid- 
IR images (see Figure [T 
G0.253 -I- 0.016 in the 



Molinari et al. (2011) place 
Elliptical ring" feature orbiting 
the Galactic center, putting it 60 — 100 pc in front of the 
Galac tic center, at a distance of ^8.4 kpc (Reid et al. 
20091. With a dust temperature of ^20 K, volume den- 
sity of >10'*cm~^, dust mass of ^10^ M©, and almost 
no signs of star-formation, GO. 253 -I- 0.016, has exactly 
the properties expe cted for the precursor to a massive. 



Arches-like cluster ([Lis et al.|[l9 94l |Lis k, Menten||1998 
Lis et al.||2001l|Ballyeral.||20l0[ ). 



Despite having the expected global properties for a 
YMC precursor, the dynamical state and fate of fu- 
ture star formation in GO. 253 -f 0.016 remains uncertain. 
While the detection of a weak water maser suggests t here 
may be a sm all number of low-mass stars forming (Lis 
et al. 1994), it is unclear whether the gas is globally 
bound and collapsing (and therefore going to form a mas- 
sive c luster) or unbound and transient (Lis & Menten 
1998 1 . This ambiguity is complicated by the extreme en- 
vironmental conditions within the Galactic center region 
(e.g. crowding, intense radiation fields, large magnetic 
fields, tidal shear). Given its importance as a potential 
precursor to a YMC, we make use of higher resolution 
and sensitivity data to study GO. 253 -I- 0.016 in detail 
and attempt to ascertain its dynamical state. 

2. SURVEY DATA AND NEW OBSERVATIONS 

2.1. Survey data 

We make use of multi-wavelength archival data 
from the United Kingdom Infrare d Digital Sky Survey 
(UKIDSS: ' ' 

Spitzer Ga 



Lawrence et al. 2007 data release 7), the 



traordinaire (GLIMPSE: [Benjamin et al |2003[ ), the Sub- 
millimetre Common-User Bolometer Array (SCUBA) 
on the James Clerk Maxwell Telescope (JCMT) ( |Di] 
Francesco et al. | 2008[) and t he methanol multi-beam sur- 
vey (MMB: [CasweU et al.||2010|). We also make use of 
ilabli 



recently available data fro m the Herschel Infrared Galac 



tic Plane Survey (HiGAL: |Molinari et al.||2"010b 
the H2O souther n Galactic plane survey (HUPS: 



et al.[2008 2011 ) and the Milli metre Astronomy Legacy 



2011) 



Walsh 



TiaS 90 GHz Survey (MALT90: | Jackson et al.|2011 



Fos- 



ter et al.|20"TT| ). Details of the surveys including observa- 



tional parameters and references outlining the data ac- 
quisition/reduction procedures can be found in Table [l] 

2.2. VLT/NACO Observations 

Deep, near-IR data were obtained with the Nasmyth 
Adaptive Optics System Near-Infrared Imager and Spec- 
trograph (NAOS-CONICA, or NACO) on the Very Large 
Telescope (VLT) on May 5th 2011. A field of view of 
1' X 1' centered on 17''42"58^ -28°42'23" (J2000) was 
imaged in J, H and Ks- Although offset from the nom- 
inal center of GO. 253 -I- 0.016, this position was chosen 
to be centered on a foreground star which was bright 
enough to use the Natural Guide Star adaptive optics 
system guiding with camera S54 and dichroic N20C80. 
After flatfielding and sky subtraction using IRAF^ stan- 
dard tools, the six dithered exposures of 10s each in each 
filter were shifted and combined into the final J, H and 
Ks images. The source extraction and PSF photometry 
were done with IRAF DAOPHOT using standard tech- 
niques and error cuts. We used Topcat^ to match all the 
stars towards G0.253 -h 0.016 observed by VLT/NACO 
with the UKIDSS data (see 



3.1) 



wavelengths were corrected to firs 



Differences in filter 
order using the prin- 
ciples describedjn Stead fc Hoare ( 2009 1 and Appendix 
A of |Sch6del et"il | ( |2010[ ). RiS5i nmg uncertainties in 
color etiects and from the statistical scatter mean the 
error in calibration is no more than ±0.1 mag at each 
wavelength. The final NACO catalog contains 358, 1001 
and 1830 sources, detected to the 5(J level up to 21.3, 
20.5 and 19.8 mag in J, H and Kg, respectively. Despite 
the error cuts, up to 2% of the sources in the catalog in 
each filter are still spurious detections. 

3. RESULTS 

Figure [1] shows continuum images of GO. 253 + 0.016 at 
3.6-8.0Aim, IQ^Ttv and 450/im. G0.253 -I- 0.016 is elon- 
gated with an aspect ratio of approximately 3:1 and is 
seen as an extinction feature from 3.6 to 70/xm with no 
obvious embedded emission sources. As such it must 
be both cold and dense and sit in front of the major- 
ity of the diffuse Galactic mid-IR background emission. 
At wavelengths >170Mm, bright emis sion is seen across 
G0.253 + 0.016 ( |Molinari et~al]|2011[ ) with similar mor- 
phology to the 450/im emission shown in Figure [l] The 
lack of embedded emission sources at wavelengths up to 
70/im combined with results from recent m aser surveys 
which show no 6.7 GHz CH3OH masers (Caswell et al. 



actic Legacy Infrared Mid-Plane Survey Ex 



^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the Nationa l Science Foundation. 
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TABLE 1 

Summary of observational survey data used in this work 



Survey /Archive 


Telcscope(s) 


A 


continuum/line 


e 


Reference (s) 


UKIDSS 


UKIRT 


1.2, 1.6, 2.2 /iin 


continuum 


<i" 


iLawrcncc ct al.l|2007| 


GLIMPSE 


Spitzer 


3.6, 4.5, 5.6, S.Oum 


continuum 


2" 


Bcniamin ct al. ( 20031 


HiGAL 


Herschel 


70, 160, 250, 350, 500 Aim 


continuum 


5-36" 


Molinari ct al. 1 2010b 201l] 


SCUBA Archive 


JCMT 


450, 850 /jm 


continuum 


8" 


(Di Francesco ct al. 20081 


MALT90 


Mopra 


3 mm 


line 


35" 


j Jackson ct al. 2011 Foster ct al.|2011 




HOPS 


Mopra 


12 mm 


line 


2' 


] Walsh ct al. 2008 2011^ 


MMB 


Parlces & ATCA 


3 cm 


line 


~1" 


(jCaswcll ct al. 201ol 



2010 



or new 22 GHz H2O maser detections (Walsh et al. 
strongly reinforces the paucity of active star tor- 



mation within the cloud. 



3.1. Near-IR extinction 

From continuum observations alone it is not possible 
to determine whether the emission from GO. 253 + 0.016 
arises from a single physical entity or from multiple 
clouds along the same line of sight. Extinction mea- 
surements can be used to distinguish between these two 
scenarios. Additionally, given the confusion and the diffi- 
culty deriving kinematic distances towards the Galactic 
center, extinction measurements also offer an indepen- 
dent distance determination. 

We used data from U KIDSS, specifically the Galac- 
tic Plane Survey (GPS: Lucas et al. 2008), and deep 
VLT/NACO observations (see § |2.2| ) to investigate the 
near-IR extinction towards GO. 253 + 0.016. We down- 
loaded the infrared sources in a 15' x 15' region cen- 
tered on GO. 253 + 0.016 from the UKIDSS data base, 
which allowed for a direct and coherent comparison of 
the cloud and off-cloud properties. The comparison of 
on- and off-source Kg vs. [H — Kg) color- magnitude di- 
agrams (CMDs) showed that there are far fewer very 
red stars in the direction of GO. 253 -I- 0.016 than in any 
of the off-source diagrams, with a clear and maintained 
deficit in the number of red stars well established by 
{H ~ Ks) = 1.7 ± 0.2mag (see top panel of Figure [2|^. 
This is equally clearly seen when comparing the NACO 
data for GO. 253 -I- 0.016 with similar field of view, al- 
beit shallower, NACO observations toward the Galactic 
Centre b y Schodcl et al (2010): Figure [2| shows the 
vs. {H - ks) CMD, and tlie histogram of (i? - X^) col- 
ors for that Galactic Centre field and for GO. 253 -I- 0.016. 
Toward the Galactic center, the red giant branch bump 
feature is clearly seen extending to {H — Kg) ~ 2.6 mag, 
after which it is effected by completeness limits. How- 
ever, despite being ^2 mag deeper, the number of sources 
toward GO. 253 -I- 0.016 drops rapidly above an {H — Kg) 
of ^1.7 mag, clearly showing that the extinction is pro- 
duced by a sharp increase in dust density as expected 
from a dense cloud, as opposed to being produced by 
several low-density clouds along the line-of-sight which 
would be seen as a much more gradual decrease in the 
number of red stars. 

We used the NACO data and Eq. 2 of [Nishiyama et ah 
( 2006 1 to estimate the distance to the cloud, using the 
red clump (RC) stars around Kg ^ 15 mag. We assume 
an absolute magnitude for the RC stars of = —1.54 



mag'*, a populati on correction AilJy = —0.07, and the 
extinction law of Schodel et al. (20101 for the Galactic 
Centre {A\ oc The distance modulus is then 

given by DM = K - Mk, + AM^f , where K is the ob- 
served, de-reddened K magnitude. We are looking for 
the distance to the near side of the cloud, so we used 
[H — Ks) = 1.7 mag to determine the extinction, since 
this is the color where we start to lose stars with re- 
spect to the Galactic Centre of Schodel et al. (2010) (see 
Figure [2j right-hand panel). Using the aforementioned 
extinction law, we derive Ak^ — 2.13 mag for an effec- 
tive wavelength of 2.168/im, which, when applied to the 
observed magnitude of the RC stars {Ks^obs = 15 ± 0.3 
mag), yields a de-reddened Ks — 12.87 ± 0.3 mag, or 
K = 12 .86 lb 0.3 considering th e difference between Ks 
and K ( [Nishiyama et al^ 2006 ) . The distance modulus 
is then DM = 14.34 ± 0.3, translating into a distance 
of 7.4 ± 1.0 kpc and independently placing it just fore- 
ground of the Galactic Centre. The greatest source of 
uncertainty in this analysis is undoubtedly the extinc- 
tion law, and we note that a difference of only 10% in 
the exponent of the extinction law in particular trans- 
lates into an approximately 10% change in the derived 
distance. 

In summary, based on the UKIDSS and NACO data we 
conclude that we have undoubtedly detected the extinc- 
tion caused by GO. 253 + 0.016, and that this is a single 
physical entity (as opposed to multiple clouds separated 
by large distances along the line of sight) with a distance 
consistent with being close in proximity to the Galactic 
Centre. Within the uncertainties, the IR-derived dis- 
tance to GO. 253 -|- 0.016 of 7.4 ± l.Okpc is consistent 
with the distance of 8.4 kpc that would be inferred if 
G0.253 -f 0.016 is part of the "1 00 pc" ring orbit i ng th e 



Galactic center, as proposed by Molinari et al.| ( |2011[ ). 
The distance to the "100 pc ring", in turn, is bas ed on 
the distance to the Galactic center determined by |Reid 
et al. (20091. As the formal uncertainty in the distance 
to the Galactic center is smaller than that from the IR- 
derived distance, we adopt a distance to GO. 253-1-0. 016 of 
8.4 kpc for the remainder of the paper. The results of the 
paper remain valid if the nearer distance were adopted. 

3.2. Dust temperature, column density, radius and dust 

mass 

HiGAL data of the inner 2° x 2° of the Galactic center 
covering the Cen tral Molecular Zone an d GO. 253 -|- 0.016 
are presented by Molinari et al. (2011). Because contin- 
uum emission at these wavelengths is typically optically- 
thin and arises from the dust, we can use the Herschel 



^ It should be noted that across the UKIDSS field examined 
there is also considerable variation in the foreground extinction, 
not surprising in a field so close to the Galactic Centre direction. 



adopt the [Schodel et al.| (I2OI0I1 
lan that of Mu = -1.59 in llMishiv 



* We 
rather than mar 01 ivi^. 
are directly comparing our data to 



value of Mf, 



-1.54 

ishiyama et alT] | |2006[ l as we 
he lormer. 
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data to derive the dust temperature and column density 
of G0.253 + 0.016. 

The dust temperature and gas column density of 
GO. 253 + 0.016 were estimated from a two-step pixel- 
to-pixel graybody fit to the Hi-GAL data. In the first 
coarse-angular-resolution step, graybody fits were calcu- 
lated using the 170 - 500 /im data, scaled to the resolution 
of the 500 ^m image (~36"). These results were used as 
inputs to the second, high-angular-resolution graybody 
fits to the 170 - 350 /im data, which were scaled to the 
resolution of the 350 fim image (~25"). Both steps used a 
fixed /? of 1.75 and did not include the 70 fim point due 
to the high optical depth at this wavele ngth. The fit- 
ting m etho d is similar to that descr ibed in [Bernard et al 



radius, R, and hnewidth, AV, through Mvir oc RAV"^, 
where the constant of proportio nality depends on the 



geometry and g as densit y profile (Ma cLaren et al 
Bertoldi & McKee 19921 IDunhani et al. 2010 



1988 



early. 



small differences in linewidth can strongly affect the 
virial mass so it is important to ensure the measured 
linewidth reflects the underlying gas kinematics and has 
not been affected via other mechanisms (e.g. different 
excitation conditions, outflows, shocks, chemistry and 
optical-depth effects). 
We used recently available data from the ongoing 



(2010) and Compiegne et al. (2011[) but with the diffuse 



comp onent of the emission removed (see Battersby et al 



MALT90 molecular line survey ( Jackson et al.|20"TT 



Fos- 



2011 for full details on the methodology] 



Figure [3] shows the derived dust temperature and col- 
umn density maps. We find the dust temperature overall 
is low, increasing smoothly from ~19K at the center to 
^27 K at the edge. There are no obvious small pockets 
of heated dust from any embedded sources. The derived 
external temperature of GO. 253 -I- 0.016 is significantly 
warmer (>35K). 

The column density and dust temperature are 
anticorrelated as expected for an externally-heated, 
dense clump. The derived peak column density is 
'-^3.3x lO^'^ cm~^ and decreases smoothly towards the 
edge. The cloud area was defined using a column density 
threshold of 3x10^^ cm~^ which covers all the emission 



ter et al. 2011 1 to investigate the velocity structure of 
GO. 253 -|- U.016. Detailed analysis and modelling of the 
emission from the many detected molecular line tran- 
sitions are underway and will be presented in a subse- 
quent paper (Rathborne et al. in preparation). Here 
we focus on the global line properties with the aim of 
deriving the most reliable, representative linewidth for 
G0.253 0.016. 

The brightest, optically-thick lines [e.g. HC0'*"(1- 
0), HNC(l-O)] show the velocity structure towards this 
region is complicated. Emission is seen from several 
distinct velocity components over the spatial extent of 
GO. 253 -t- 0.016 defined in § |3.2| The most prominent of 
these are: a component at '-^O kms~^ seen at the north- 
east edge of the filament; a component at ~35kms~^ 
covering a similar emission area as the spatial extent 
of GO. 253 -I- 0.016 defined above; and a component at 
^75kms^^ seen towards the south west edge of the fila- 



is 



at T< 27 K. Within this area the average column density ment. Given the spatial extent of the 35kms~ compo- 
nent closely matches that of the dust emission, we take 
this component as the emission from GO. 253 + 0.016, 
with the other components being from unrelated clouds 
along the line of sight. 

The 35kms~^ component shows evidence for shifts in 
the peak Visr with changes in spatial position (ie a veloc- 
ity gradient). This means that taking a spectra from an 
individual spatial position would not accurately recover 
the integrated linewidth over the whole region, as the 
linewidth-broadening from the velocity gradient would 
not be taken into account. To overcome this requires 
extracting spectra which are integrated over the spatial 
extent of G0.253 + 0.016. 

We selected bright, optically-thin emission as the most 
likely reliable tracer of the underlying gas kinematics. 
The top panel of Figure |4] shows Hanning-smoothed spec- 
tra from the H"CO+(1-0) and HN"C(l-0) transitions 
integrated over the spatial extent of the cloud (as de- 
termined in 



1x10 cm~ . The emission is elongated with semi- 
major and semi-minor axes of l.'9x0.'7, corresponding to 
a physical radius of 4.7xl.7pc for a distance of 8.4 kpc. 
Taking the geometric mean, we derive an effective radius 
of 2.8pc for GO. 253 -I- 0.016. 

Given the strong dependence of dust emissivity with 
temperature and to account for the dust temperature 
variation across GO. 253 -|- 0.016, we calculate the dust 
mass directly from the derived column density rather 
than from mm continuum emission which assumes a 
single dust temperature. The dust mass was calcu- 
lated by multiplying the column density of each pixel 
by the pixel area, assuming a mean mo lecular weight of 
2.8xMH-atom (Kauffmami et al. 2008) and gas-to-dust 
ratio of 100:1 ("Hildebrand 1983^ and summing over the 
pixels in the cloud area. The total mass is not particu- 
larly sensitive to the column density threshold, varying 
by <10% for threshold values of 1— 5x 10^^ cm~^. We de- 
rive a dust mass for G0.253 + 0.016 of 1.3xlO^M0 and 
estimate the uncertainty to be of order 20%. 



3.2). The line profiles for both transi- 



3.3. Molecular line data: linewidth and virial mass 

To determine whether GO. 253 + 0.016 is 
gravitationally-bound we compared the dust mass with 
the virial mass, Mv 



^ Due to the metallicity gradient in the disk jBalser et al.|201l| |, 
it is possible the gas-to-dust ratio in the CMZ may be lower than 
the canonical value of 100:1 commonly adopted in the literature. 
However, as individual metallicity measurements across t he CMZ 
vary from 1 to ~4 X Solar (Shields & Fcrland 1994, Macd a et al.| 
[2002|[Naj arro et al.|2009^ and no independent measurements exist 
tor ciu.253 -I- O.UIB, we opt for the standard gas-to-dust ratio of 
100:1, while noting the dust-based mass may possibly be lower. 



tions are similar, giving confidence that they are ro- 
bustly tracing the underlying gas kinematics. Both the 
Visr = 35kms~^ component from GO. 253 -I- 0.016 and the 
Visr — 0kms~^ component are seen but the 75kms~^ 
component is too faint in these optically-thin transitions 
to be detected. The bottom panel of Figure |4] shows the 
spatially-averaged HN^'^C(l-O) emission with no Hanning 
smoothing. This emission was fit with a two-component 
Gaussian profile using the tasks within the CLASS^ soft- 
ware package. To avoid biasing the fit results we did not 
constrain any of the parameters but rather left them all 

^ CLASS is part of the GIL DAS software environment: 
|http://iram.fr/IRAMFR/GILDAS/| 
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as free variables. From the resulting fits, which are over- 
layed on the spectra in the bottom panel of Figure [4j we 
derive the peak velocity and FWHM of the gas associ- 
ated with G0.253 -I- 0.016 to be V^r = 36.1 ± 0.4kms-i 
and AV= 15.1 ± l.Okms"^, respectively. We adopt an 
upper limit linewidth of 16kms~^ for GO. 253 + 0.016 for 
the remainder of the paper. 

The derived virial mass can vary by an order of mag- 
nitude depending on the assumed geometry and density 
distribution. Assuming a spherical source with a den- 
sity distrib ution oc r gives a v irial mass estimate of 
9x10'*Mq ( [MacLaren et al.|l988 [ Eq 3). However, more 
realistically, for a centrally-condensed cloud elongated 
with axis ratio 3:1 a nd densitv oc r~^'^, the v i rial mass 
would_be4xlO^M0 ( iBertoldi fc McKee||l992| [Dunham" 
et al.|2010 ) . Despite the assumptions and systematic un- 
certainties, the dust and virial mass estimates agree to 
within a factor of a few. Thus it is likely GO. 253 -|- 0.016 
is close to virial equilibrium. 

4. DISCUSSION 

The global properties of GO. 253 -I- 0.016 derived in § |3] 
are summarized in Table [2] Based on these results, 
GO. 253 -I- 0.016 is a strong candidate precursor of an 
Arches-like YMC. The only other two similar Galactic 
YMC-forming clouds that have been studied to date are 



Sgr B2 and W49 A (e.g. Goldsmith et al. 1990 Alyes 



fc Homeier||2003|) . However, these are both much more 



evolved, with massive cluster formation well underway. 
The powerful feedback from ongoing star formation has 
strongly affected the cloud structure in these regions, so 
observations of the gas do not probe the initial conditions 
prior to the onset of star formation. 

We speculate that GO. 253 -I- 0.016 may be unique in 
representing the initial conditions of a Galactic molecu- 
lar cloud on the verge of forming a YMC. As such, its 
detailed study can reveal important clues about massive 
cluster formation and help test theoretical models. Cur- 
rent theories suggest that the early evolution of stellar 
clusters depends cruciall y on the spatial and kinematic 
distribution of stars (e.g. Goodwin & Bastian 2006; Alli- 
son et al.|2010 1 . However, while the distribution of emcr- 
gent stellar populations as "hierarchical" or "centrally- 
condensed" may be straightforward to quantify when an- 
alyzing stars, quantifying the small-scale distrib ution of 



the gas in a similar way is complicated (see e.g. Lomax 
et al. 20111. Because gas can not easily be counted in 
discreet, observationally-defined units, the results of al- 
gorithms attempting to quantify gas substructure are no- 
toriously subjectiye_to_t^ of user-defined inputs 
(see e.g. |Pineda et ar]|2009[ ). 

Regardless of the exact details of the formation mech- 
anism, the properties of the emergent stellar population 
must be directly related to the initial global gas prop- 
erties. Assuming a reasonable star f ormation efficiency 
(SFE) (e.g. 30% for gas > 10^ cm'^; [Xlves et al.|[2007| , 
GO. 253 + 0.016 may form a cluster with stellar mass and 
density similar to that of other YMCs like the Arches. 
Whether or not such a cluster remains gravitationally- 
bound depends on the SFE and environmental condi- 



tions (tidal shear, interaction with GM Cs etc, e.g. Krui 
jssen 2011 Elmegreen & Hunter 2010). However, it is 
clear that GO. 253 -I- U.016 has enough mass to form an 
Arches-like cluster without accreting any additional gas 



or stars from outside the present-day observed boundary. 

With only a single candidate YMC precursor it is not 
possible to determine how representative GO. 253 -|- 0.016 
is of all YMC progenitors. It is clearly desirable to search 
for other proto-YMC molecular clouds, especially within 
the Galaxy. In the last decade much observational effort 
has been directed towards searching for, and characteriz- 
ing the physical properties of molecular clouds which are 
likely the progenitors of stellar clusters. As a result, a 
large fraction of the gas in the Galaxy available to form 
YMCs has been surveyed in dust continuum and dense 
molecular gas tracers. 

Figure [5] summarizes the results from two of these sur- 
veys and places GO. 253 -|- 0.016 in the context of Galac- 
tic molecular clouds and stellar clusters. Plus symbols 
show ammonia clumps detected in HOPS (Walsh et al. 
2011 1. The critical density of ammonia is close to th e vol- 
ume density threshold of ^10'* cm~'^ proposed by Lada 



et al.] (2010) as the threshold required to form stars. As 
such HUPS is ideal for identifying and characterizing the 
properties of star forming gas within the Galaxy. Green 
crosses show IRDCs from the survey of [Rathborne et aL] 
(|2006j). IRDCs are thought to be precursors of high mass 
stars and clusters. The region of the mass-radius plane 
covered by the HOPS clouds and IRDCs is similar to 
that of clouds dete cted in other surveys, bot h targeted 
(e.g. as compiled by'KaufFmann fc Pillai|20l0 | and blind, 
large-area s urve ys (BGPS - Aguirre et al. "(2011); AT- 



LASGAL - [Schuller et al.| ( |2009[ ), I'achenberg et al sub 



mitted and Beuther et ah priv. comm). The hatched 
rectangles s how the mass-radius range o f different stel- 



lar clusters (Portegies Zwart et al. 2010) and the black 
dots show Galactic YMCs. With tne exception of a few 
clouds which may form small YMCs, assuming a rea- 
sonable star formation efficiency, most of the observed 
molecular clouds seem destined to form open clusters. 
GO. 253 + 0.016 is marked with a red star and clearly 
stands out. It is significantly more massive for its size 
compared to the other Galactic molecular clouds. It is 
potentially the only object with an observed gas mass 
and density which may be able to directly form a clus- 
ter with stellar mass and densities in the YMC/globular 
cluster regime. 

The potential uniqueness of GO. 253 + 0.016 has pro- 
found implications for YMC formation in the Galaxy. We 
posit three scenarios. If all YMCs form from molecular 
clouds with similar properties to GO. 253 + 0.016, either 
i) GO. 253 -I- 0.016 is unique in the Galaxy, or, ii) surveys 
have missed the other similar clouds. Alternatively, iii) 
if the properties of GO. 253 -I- 0.016 are not representa- 
tive of other molecular clouds destined to form YMCs, 
there may be precursors to YMCs in the Galaxy which 
have already been observed but not identified as YMC 
precursors. 

To distinguish between these scenarios it is important 
to predict how many proto-YMC molecular clouds are 
likely to exist at any given time in the Galaxy. One 
w ay to es t imate this is following the analytical approach 
of Gieles (2009). Assuming a star formation rate of 1 
M0/yr with 10% of stars forming in clusters and a clus- 
ter mass function index of —2, we would expect ~10 
clusters of >1O''M0 with ages <10Myr. This estimate 
is consistent with the number of YMCs currently known 
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in the Galaxy. If the number of Galactic YMCs is con- 
stant, assuming it takes of order a Myr to evolve from 
precursor to forming cluster we would expect from zero 
to a few precursors at any given time. An alternative 
method to estimate the expected number of proto-YMC 
molecular clouds in the Galaxy is extrapolating the ob- 
served embedded cluster mass function (ECMF), which, 
for embedded clusters up to 10'^ M f^ within 2 kpc of the 
Sun, is a power- law with index —2 (Lada & Lada 2003). 
Assuming this represents the Galactic ECM.b' (see.GieleS| 
[2009^ , one can scale it to the total volume of the Galac- 
tic disk to estimate the number of embedded clusters ex- 
pected in the Galaxy. Extrapolating this to higher clus- 
ter masses, we expect ~ 25 young clusters with masses 
between 10^ and 10^ M©. Assuming a cluster lifetime 
of 10 Myr we predict ~ 2.5 proto-YMCs per Myr in the 
Galaxy. Both of these approaches suggest there may be 
other Galactic YMC progenitors. While these estimates 
are approximate, they predict that only a small num- 
ber of proto-YMCs should exist in the Galaxy at any 
given time. However, the uncertainties are sufficiently 
large that it is not possible to use these estimates to 
rule out GO. 253 + 0.016 as being unique. Therefore, it is 
currently not possible to distinguish between scenarios i) 
and ii). 

Returning to scenario iii), there are many molecu- 
lar cloud complexes known with sufficient mass to form 
YMCs but their present-day gas density is much lower 
than the expected final YMC stellar density. Examples 
of these are seen clearly in Figure [5] as clouds with mass 
> 10^ Mq and radii >^ 8 pc. Such clouds seem destined 
to form massive, but unbound, associations. However, 
it may be possible to form the required YMC stellar 
densities over time from the convergence of molecular 
gas flows on large scales and subsequent global gravi- 
tational collapse. Indeed, such large-scale gravitational 
collapse has been directly observed towards a number of 
mass ive proto clusters (e.g. W33A dCalvan-Madrid et al. 



2010 



DR21 ([Schne ider et itlJ^mr GTUWl^eto et al. 
Baobab Liu et al. 20^7"'^67 ( [Longmore et al. 



20111) and is predicted by numerical si mulations (Smith 



erar|2009[[Klessen fc Hennebelle|20To[ ) . In this scenario, 



star formation would proceed over several crossing times 
leading to a large age spread in cluster members. In this 
regard, it is interesting to note that an age spread of 
> 10 Myr among cluster members are reported in sev- 
eral YMCs (NGC 3603, NGC 346 30 Doradus - [Bee 



Ol [De Marchi et al.[[2011a|b| ). However 
aebate in the literature over the uncer- 



cari et al.[[201 0| 

there is some 

tainty in the age determination of cluster members using 
this method and whether or not thes e data are also con- 
sistent with age sp r eads of <3M yr (iTobin et al. 2009[ 
Baraffe et all 20091 ILittlefair et al. 2011: Jettrics et al 



2011 r Detailed studies of the large-scale gas motions 
in the youngest and most massive molecular clouds are 
needed to determine whether the gas/stellar density will 
increase over time to form YMCs as opposed to more 
diffuse and gravitationally-unbound OB associations as 
would be predicted by their current global gas density. 

With this scenario in mind, it is interesting to note that 
large-scale emission from shocked-gas tracers is detected 
towards GO. 2 53 -I- 0.016. Co mbined with gas tempera- 
tures of 80 K (Lis et al.[[2001 Rathborne et al. in prep). 



which is much warmer than the dust temperature, this 
suggests GO. 253 -I- 0.016 may have formed from a cloud- 
cloud collision or convergence of large scale flows. As the 
gas has not yet had the chance to heat up the dust, this 
must have happened recently. If this is the case, before 
the postulated collision/convergence, the gas properties 
of GO. 253 + 0.016 may have been much more similar to 
the typical Galactic molecular cloud population and thus 
may not have stood out as extreme in Figure [5[ 

Given the location of GO. 253 -I- 0.016 near to the Galac- 
tic center, it is tempting to invoke the extreme con- 
ditions at the Galactic center (e.g. interstellar radia- 
tion fleld enhanced by 10"^, external pressure P/k~10*, 
strong magnetic fields, gamma rays) to justify why it 
may be unique. In this regard it is interesting to note 
that GO. 253 + 0.016 lies at a projected distance of only 
18 pc from the Arches and Quintuplet clusters. While 
it is not possible to directly link the formation of the 
Arches with GO. 253 -I- 0.016, it seems plausible that the 
Arches could have formed at its present distance from 
the Galactic center. On larger scales GO. 253 -I- 0.016 ap- 
pears to be part of a 10^ Mq filament of interconne cted 
clumps of w hich GO. 253 -I- 0.016 is the most massive (Lis 
et al.[[200T ). This filament is itself thoug ht to be part of 
the 3x10' Mq, 100 pc ring identified by Mohnari et al. 



(2011 1. The 100 pc ring is orbiting the Galactic center at 
^BOkms"^ so a single orbit would take ^3.6 Myr. Given 
the free-fall time for G0.253 + 0.016 is <lMyr (see Ta- 
ble [2|, it seems unlikely GO. 253 -I- 0.016 would survive a 
whole orbit without additional support - not implausible 
given the extreme conditions at the Galactic center. If 
such additional support could be maintained against the 
rapid turbulent energetic decay, this may explain why 
the cloud has remained dense and starless. 

Alternatively, the fate of the GO. 253 -I- 0.016 may be 
more intimately linked to the dynamics of the 100 pc 
ring. Gas is thought to enter the ring at the intersec- 
tion points of the innermost stable XI orbits and the X2 
orbits. These intersection or "crashing" points coincide 
with Sgr B2 and Sgr C - the only places in the 100 pc ring 
with prodigious star formation activity. If the gas com- 
prising GO. 253 -|- 0.016 were to have entered the 100 pc 
ring at Sgr C it would have joined only a small fraction of 
an orbit time ago. In around a free-fall time, the whole 
10^ Mq filament in which GO. 253 -I- 0.016 is located will 
have reached the position of Sgr B2. As such, this larger 
filament could be in the process of forming a 10^ Mq, 
Sgr B2-like complex. 

However, while it is tempting to invoke the ex- 
treme conditions at the Galactic center to justify why 
GO. 253 + 0.016 may be unique, several of the most mas- 
sive YMCs in the Galaxy are not located near the Galac- 
tic center suggesting the extreme environmental condi- 
tions are not a necessary condition to form a YMC. 

5. CONCLUSIONS 

We make use of existing multi-wavelength data includ- 
ing recently available far-IR continuum (Herschel/Hi- 
GAL) and mm spectral fine (HOPS and MALT90) 
data and present new, deep, multiple-filter, near-IR 
(VLT/NACO) observations to study the infrared-dark 
cloud GO. 253 -I- 0.016. From these data we have de- 
rived the global properties of GO. 253 -I- 0.016. It has 
a high mass (Mdust ~ 1-3 x 10^ Mq), low dust tern- 
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perature (~23K), small radius ('^2.8pc), high vol- 
ume density (7.3x 10^ cm""^) and high column density 
(3.5xl0^^cm~^). It is close to virial equilibrium with 
almost no signs of active star formation. As such it ap- 
pears to be a prime candidate for the initial conditions of 
a molecular cloud destined to form an Arches-like YMC. 

Comparing the properties of GO. 253 + 0.016 to other 
Galactic dense molecular clouds shows it to be extreme. 
We discuss implications of this for the formation of mas- 
sive protoclusters and posit three scenarios. If all YMCs 
form from molecular clouds with similar properties to 
G0.253 + 0.016, cither i) G0.253 + 0.016 is unique in the 
Galaxy, or, ii) surveys have missed the other similar 
clouds. Currently it is not possible to distinguish be- 
tween scenarios i) and ii). Nevertheless, clouds with 
properties like GO. 253 -I- 0.016 must be very rare in the 
Galaxy making GO. 253 + 0.016 extremely important for 
testing massive cluster formation models. Alternatively, 
in scenario iii) the properties of GO. 253 + 0.016 do not 
represent those of other molecular clouds destined to 
form YMCs. We note there are many molecular clouds in 
the Galaxy with sufficient mass to form YMCs but their 
present-day gas density is much lower than the expected 
YMC stellar densities. However, it may be possible to 
form the required stellar densities over time from the con- 
vergence of molecular gas flows on large scales and subse- 
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quent global gravitational collapse, as observed towards 
several lower-mass protoclusters. Detailed studies of the 
large-scale gas motions in the youngest and most mas- 
sive molecular clouds arc needed to determine whether 
the gas/stellar density will increase over time to form 
YMCs as opposed to more diffuse and gravitationally- 
unbound OB associations as would be predicted by their 
current global gas density. 

In future work we will study the internal structure 
and kinematics of GO. 253 -|- 0.016 in detail to investigate 
its recent dynamical history, test the plausibility that it 
formed from cloiid-cloud collisions or converging large- 
scale flows, and directly test models of massive proto- 
cluster formation. 
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LINEWIDTH (AV), volume DENSITY (p), COLUMN DENSITY (N^f^), CLOUD-CROSSING TIME (tec), SOUND-CROSSING TIME (tsc) AND FREE-FALL 

TIME (tff). 



M 


D 


R 


Tdust 


AV p 


Nff2 tec 


^SC 




(Mq) 


(kpc) 


(pc) 


(K) 


(km/s) (cm^'*) 


(cm"2) (Myr) 


(Myr) 


(Myr) 


1.3E5 


8.4 


2.8 


19-27 


16 7.3E4 


3.5E23 0.17 


8 


0.74 




Littlefair, S. P., Naylor, T., Mayne, N. J., Saunders, E., & 

Jeffries, R. D. 2011, MNRAS, 413, L56 
Lomax, O., Whitworth, A. P., & Cartwright, A. 2011, MNRAS, 

412, 627 

Longmore, S. N., Burton, M. G., Barnes, P. J., & et al. 2007, 

MNRAS, 379, 535 
Longmore, S. N., Pillai, T., Keto, E., Zhang, Q., & Qiu, K. 2011, 

ApJ, 726, 97 

Lucas, P. W., Hoare, M. G., Longmore, A., & et al. 2008, 

MNRAS, 391, 136 
MacLaren, I., Richardson, K. M., & Wolfendale, A. W. 1988, 

ApJ, 333, 821 

Maeda, Y., Baganoff, F. K., Feigelson, E. D., Morris, M., Bautz, 
M. W., Brandt, W. N., Burrows, D. N., Doty, J. P., Garmire, 
G. P., Pravdo, S. H., Ricker, G. R., & Townsley, L. K. 2002, 
ApJ, 570, 671 

Molinari, S., Bally, J., Noriega-Crespo, A., & et al. 2011, ApJL, 
735, L33-I- 

Molinari, S., Swinyard, B., Bally, J., & et al. 2010a, A&A, 518, 
L 100-1- 

— . 2010b, PASP, 122, 314 

Najarro, F., Figer, D. F., Hillier, D. J., Geballe, T. R., & 

Kudritzki, R. P. 2009, ApJ, 691, 1816 
Nishiyama, S., Nagata, T., Sato, S., & et al. 2006, ApJ, 647, 1093 
Peretto, N. & Fuller, G. A. 2009, A&A, 505, 405 
Pillai, T., Wyrowski, F., Carey, S. J., & Menten, K. M. 2006, 

A&A, 450, 569 

Pineda, J. E., Rosolowsky, E. W., & Goodman, A. A. 2009, 

ApJL, 699, L134 
Portegies Zwart, S. F., McMillan, S. L. W., & Gieles, M. 2010, 

ARA&A, 48, 431 
Purcell, C. R., Balasubramanyam, R., Burton, M. G., & et al. 

2006, MNRAS, 367, 553 
Purcell, C. R., Longmore, S. N., Burton, M. G., & et al. 2009, 

MNRAS, 394, 323 
Rathborne, J. M., Jackson, J. M., & Simon, R. 2006, ApJ, 641, 

389 



Rathborne, J. M., Johnson, A. M., Jackson, J. M., Shah, R. Y., & 

Simon, R. 2009, ApJS, 182, 131 
Reid, M. J., Menten, K. M., Zheng, X. W., & et al. 2009, ApJ, 

700, 137 

Santangelo, G., Testi, L., Gregorini, L., & et al. 2009, A&A, 501, 
495 

Schneider, N., Csengeri, T., Bontemps, S., & et al. 2010, A&A, 
520, A49+ 

Schodel, R., Najarro, F., Muzic, K., & Eckart, A. 2010, A&A, 
511, A18+ 

SchuUer, F., Menten, K. M., Contreras, Y., & et al. 2009, A&A, 
504, 415 

Shields, J. C. & Ferland, G. J. 1994, ApJ, 430, 236 

Simon, R., Jackson, J. M., Rathborne, J. M., & Chambers, E. T. 

2006a, ApJ, 639, 227 
Simon, R., Rathborne, J. M., Shah, R. Y., Jackson, J. M., & 

Chambers, E. T. 2006b, ApJ, 653, 1325 
Smith, R. J., Longmore, S., & Bonnell, L 2009, MNRAS, 400, 

1775 

Sridharan, T. K., Beuther, H., Schilke, P., Menten, K. M., & 

Wyrowski, F. 2002, ApJ, 566, 931 
Stead, J. J. & Hoare, M. G. 2009, MNRAS, 400, 731 
Swinbank, A. M., Small, L, Longmore, S., & et al. 2010, Nature, 

464, 733 

Taylor, M. B. 2005, in Astronomical Society of the Pacific 
Conference Series, Vol. 347, Astronomical Data Analysis 
Software and Systems XIV, ed. P. Shopbell, M. Britton, & 
R. Ebert, 29—1- 

Tobin, J. J., Hartmann, L., Furesz, G., Mateo, M., & Megeath, 

S. T. 2009, ApJ, 697, 1103 
Walsh, A. J., Breen, S. L., Britton, T., & et al. 2011, MNRAS, 

416, 1764 

Walsh, A. J., Lo, N., Burton, M. G., & et al. 2008, PASA, 25, 105 
Whitmore, B. C. 2000, ArXiv Astrophysics e-prints 



GO. 253 + 0.016: a molecular cloud progenitor of an Arches-like cluster 



9 



11.0- 
11.5 - 
12.0 
12.5 ■ 
13.0 - 

9 

" 13.5 - 
14.0 ■ 
14.5 ■ 
15.0 - 

15.5 - 

16.0 
0.0 

























— 








• 






• 












































— ■ 


- • 




: 


•a ■ 








— - 






'i 


, * 
















• * 











































11.0 
11.5 
12.0 
12.5 
13.0 
' 13.5 
14.0 
14.5 
15.0 
15.5 



2.0 
(H-Ks) 



14 



Ml :.' 



































-* — 












• 












• • 






























. * *•* 








































































0.5 1.0 1.5 2.0 2.5 3.0 3.5 4. 

(H-Ks) 



1000 




1 2 

(H-Ks) 

Fig. 2. — [Top panel] Kg vs. {H — Ks) color-magnitude diagrams for UKIDSS sources within a circle of radius 45" centered on 
GO. 253 + 0.016 [left] and an offset position at I = 0.32°, 6 = -0.02° [right]. Although the UKIDSS data are confusion-limited 
in sensitivity and not as deep as the NACO data (bottom panel) a clear decrease in source counts towards GO. 253 -I- 0.016 is 
seen at H — Ks ~ 1.7 ± 0.2 mag. [Bottom left panel] Ks vs. ( H — Ks) c olor- magnitude diagram. The black dots show sources 
detected in the VLT/NACO data towards GO. 253 -I- 0.016 (see §|2^&|3l]). The yellow dots and grey s cale sh ow the VLT/NACO 



sources and number density, respectively, detected towards the Galactic center from |Schodel et al.| ( |2010| . Both VLT/NACO 
observations have the same 1' x 1' field of view but the GO. 253 -I- 0.016 obs ervations are ~2mag ~aeeper (as illustrated by the 



dotted and dashed lines showing the approximate 5a sensitivity limits of the Schodel et al. (20101 and GO. 253 -f- 0.016 datasets 
respectively). In the Schodel ct al. (2010J data, the red giant branch bump teature continues to {H — Ks) ~ 2.4 mag and is 
still visible at 2.6 mag, where it becomes strongly effected by completeness limits. Despite being ~2mag deeper and looking 
along similar lines of sight (and hence similar foreground extinction and expected stellar populations), the red giant branch 
bump fea ture is sharply trun cated in the GO. 253 -I- 0.016 observations. [Bottom right panel] Histogram of the ( H — Ks) colors 
for the Sc hodel et"aL ] (j2010| catalog [black], and for the NACO sources in GO. 253 + 0.016 brighter than the [Sch odel et al. 
([2010) detection limit [rcdj. Both distributions are similar for {H — Ks) < 1.7 but the number of GO. 253 -I- 0.016 sources shows 
a sharp cutoff at redder colors, despite the observations being ~2 mag deeper. This is clear evidence that we have detected 
the extinction from GO. 253 -I- 0.016, and that this is a single entity as opposed to multiple clouds along the line of sigh t. T he 
distance determined from these data is consistent with GO. 253 -I- 0.016 being in close proximity to the Galactic center (§ |3.1[ ). 
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Fig. 3. — Dust temperature [left] and column density [right] maps of GO. 253 + 0.016 derived from the Hi-GAL data. The images are 
12pc on a side for the distance of 8.4kpc. The temperature contours are 19, 20, 28K and the column density contours are 5, 10, 
35 X lO^'^ cm~'^. T he c olor scale shows all pixels above the threshold column density of 3xl0^^cm~'^ used to define the boundary of 



GO. 253 + 0.016 (see § |3.2[ |. The dust temperature overall is low, increasing smoothly from ~19K at the center to ~27K at the edge. 
There are no obvious small pockets of heated dust from any embedded sources. The derived external temperature of GO. 253 + 0.016 is 
significantly warmer (>35K). The column density and dust temperature are anticorrelated as expected for an externally- heated, dense 
clump. The derived peak column density is ~3.3xl0'^'^ cm~^ which decreases smoothly towards the edge. The average column density in 



the region above the threshold cutoff is ~lxlO^ 
mass to be 2.8pc and 1.3x10^ Mq, respectively. 



cm . Based on the column density maps we derive GO. 253 + 0.016's effective radius and 
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Fig. 4.— [Top] Hanning-smoothed H13cO+(1-0) [blue] and HNi^C(l-O) [red] spatially -averaged spectra across GO. 253 + 0.016 from 
the MALT90 survey. Emission from both these transitions is expected to be optically-thin so should reliably trace the underlying gas 
kinematics. The line profiles for both transitions are similar, with two components: a brighter one at Visr = 35kms~^ from GO. 253 -I- 0.016 
and a weaker one at Via^ = Okms~^ from an unrelated cloud along the line of sight. [Bottom] The black line shows the same HN^^C(l-O) 
line profile as above, but without any Hanning smoothing. The green line shows the results of a two-component Gaussian fit to this profile, 
from which a line FWHM of 15.1 it l.Okms" is derived for the gas associated with GO. 253 -I- 0.016. 
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Fig. 5,— Radius versus mass for Galactic dense, cluster-forming molecular clouds. Plus symbols show ammonia clouds detected in 
HOPS (Walsh et al. 2011} (blue/ brown denote an assume d near/far kinematic distance, respectively). Green crosses show infrared dark 
clouds (IRDGs) fr om the survey of |Rathborne et al.|(|2006| . The hatched rectangles show the mass-radius range of different stellar clusters 
portegies ZwariTc t al. 2010). The black dots siiow Galactic young massive clusters. With the exception of a few clouds which may form 
small YMGs, assuming a reasonable star formation efficiency, most of the observed molecular clouds seem destined to form open clusters. 
GO. 253 + 0.016 is marked with a red star and clearly stands out as unique. It has a mass and radius that would be expected of a molecular 
cloud progenitor of a large YMC or a globular cluster. The dashed lines show constant density and free-fall time. 



